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Abstract 

Exact chiral symmetry at finite lattice spacing would preclude the axial anomaly. In order to 
describe a continuum quantum field theory of Dirac fermions, lattice actions with purported exact 
chiral symmetry must break the flavor-singlet axial symmetry. We demonstrate that this is indeed 
the case by using a minimally doubled fermion action. For simplicity we consider the Abelian axial 
anomaly in two dimensions. At finite lattice spacing and with gauge interactions, the axial anomaly 
arises from non-conservation of the flavor-singlet current. Similar non-conservation also leads to 
the axial anomaly in the case of the naive lattice action. For minimally doubled actions, however, 
fine tuning of the action and axial current is necessary to arrive at the anomaly. Conservation of 
the flavor non-singlet vector current additionally requires the current to be fine tuned. Finally we 
determine that the chiral projection of a minimally doubled fermion action can be used to arrive at 
a lattice theory with an undoubled Dirac fermion possessing the correct anomaly in the continuum 
limit. 
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I. INTRODUCTION 



The study of strongly interacting gauge theories has made remarkable progress due to 
numerical computations on Euclidean space-time lattices. Lattice gauge theory provides a 
gauge invariant regularization of vector-like theories, such as QCD. Inchision of a desired 
number of light fermionic degrees of freedom, however, is not a simple task. The naive dis- 
cretization, D, of the massless free fermion Dirac operator maintains exact chiral symmetry, 
{£)~^,75} = 0, but produces fermion doubling. For each naive fermion in (i-dimensions, 2*^ 
Dirac fermions emerge in the continuum limit. Under some reasonable assumptions, one 
can make this result more general. The Niclscn-Ninomiya theorem [1] essentially states 
that chiral symmetry must be sacrificed to avoid doubling. This was the approach taken 
in the pioneering work of Wilson [2]. The Wilson Dirac operator, Dw, avoids fermion dou- 
bling because it contains an irrelevant operator that breaks chiral symmetry, so that one 
has {1)^^,75} 7^ 0. At some level, chiral symmetry must be broken by all lattice actions 
with any hope of describing Dirac fermions in the continuum limit. The reason is that exact 
chiral symmetry precludes the flavor-singlet axial anomaly. For the Wilson action, chiral 
symmetry breaking introduced by the Wilson term leads to the correct axial anomaly in the 
continuum limit [3] (we provide a simple derivation in the Appendix). 

Rather elegant solutions to the fermion doubling problem exist. Starting from continuum 
QCD and smearing the fermion fields around lattice sites, one can deduce the properties 
of an ideal lattice Dirac operator, Dqw- Such an operator satisfies what is known as the 
Ginsparg- Wilson relation [4], namely {Dq]^,'^^^ — 75. While Ginsparg- Wilson fermions 
do not have the usual continuum chiral symmetry, they maintain a lattice version of chiral 
symmetry, which guarantees good chiral symmetry properties of the theory; but, at the same 
time, ensures that the correct axial anomaly is recovered [5]. Domain wall fermions [6], and 
their four-dimensional reduction, overlap fermions [7], are explicit solutions to the Ginsparg- 
Wilson relation.^ While elegant, these fermions are numerically costly. It is of great interest 
to have inexpensive chiral fermions for QCD computations, and additionally for studies of 
strongly coupled dynamics in theories beyond the standard model, for a review of such 
theories, see [9]. 

The naive action has a form of chiral symmetry not encompassed by the Ginsparg- Wilson 
relation. As the low-energy modes emerge from several regions within the Brillouin zone, 
the fiavors of naive fermions (often referred to as tastes) are non-local in space-time.^ While 
the free naive action maintains an exact fiavor-singlet axial symmetry, the corresponding 
axial transformation of the fermion field is non-local. In the presence of gauge interactions, 
the flavor-singlet axial symmetry must be broken. Consequently the correct axial anomaly 
is recovered in the continuum, as was realized in the early days of lattice QCD [10]. Today 
one commonly encounters the statement that naive fermions do not have an axial anomaly, 
see, for example, [11]. Such statements are solely based on the local axial transformation, 

— >■ exp(i^75)'03., which is a singlet transformation of the naive fleld. This transformation 
generates an exact symmetry of the interacting theory; but, in terms of a flavor interpreta- 
tion, careful consideration reveals that it is non-singlet [3, 10]. 

Here we detail the way in which the axial anomaly emerges for theories with chiral 

^ The Ginsparg- Wilson relation can also be realized through the construction of perfect lattice actions [8] . 

^ While there is no unique way to define the flavors of naive fermions, all choices are non-local and differ 
by lattice-spacing artifacts. The mildest non-locality possible has a range of just one lattice spacing. For 
this choice, the relative coordinates within an elementary hypercube can be regarded as internal degrees 
of freedom which give rise to the flavors. Rewritten in this way, the action becomes local and contains 
flavor-breaking interactions. q 



symmetry of the form {D~^,75} = 0. In our study, we use the minimally doubled fermion 
action of Wilczek [12].^ This simplifies the flavor content of the continuum theory compared 
to the naive fermion action, and consequently simplifies the discussion of flavor singlet and 
non-singlet currents. Recently interest has been renewed in minimally doubled fermion 
actions due to graphene-inspircd constructions [14]. These actions, however, maintain fewer 
symmetries over the Wilczek action [15], giving us another reason to focus on it. As with the 
naive action, the local axial transformation, ip^ exp{i9j^)ipx, generates a conserved flavor 
non-singlet current. On the other hand, the flavor-singlet axial transformation of the fermion 
fleld is inherently non-local (symmetry transformations for minimally doubled fermions are 
collected in Table I). In the presence of gauge interactions, the corresponding singlet current 
is not conserved, and is the source of the chiral anomaly. An important additional finding of 
our study is that the non-singlet vector current is not conserved. Current operator must be 
fine tuned in order to recover conservation of the non-singlet vector current in the continuum 
limit. A similar malady affects the divergence of the singlet axial current; hence, fine tuning 
is also required to arrive at the chiral anomaly. We also analyze the chiral projection of the 
Wilczek action, and demonstrate how a single Dirac fermion emerges from its continuum 
limit. 

We organize our presentation as follows. The free minimally doubled action is spelled 
out in Sec. II, along with its symmetries, both discrete and continuous, and the form of 
the conserved vector and axial-vector currents. The inclusion of gauge interactions is in- 
vestigated in Sec. III. When gauge interactions are included, the flavor-singlet axial-vector 
current is no longer conserved. We use this observation to calculate the chiral anomaly in 
Sec. IV. where, for simplicity, we consider only the two-dimensional case. The flavor non- 
singlct vector symmetry is also violated by gauge interactions. We show that conservation 
of the corresponding current requires fine tuning of a marginal operator. There is a similar 
such marginal operator for the fiavor-singlet axial current, and it must be tuned to arrive 
at the correct chiral anomaly. Chiral projection is taken up in Sec. V, where we consider 
the minimally doubled action for two-component Weyl fermions. In the continuum limit, 
we show how this theory can describe an undoubled Dirac fermion with the correct chiral 
anomaly. The Appendix gives a simple derivation of the chiral anomaly for the case of the 
Wilson lattice action. Finally, we summarize our findings in Sec. VI. 



II. MINIMALLY DOUBLED FERMIONS 
A. Free action and its symmetries 

To begin our discussion, we consider a free lattice theory of minimally doubled fermions. 
In the absence of gauge interactions, the Wilczek action is given by 

X,IJ. x,j 

where sums over /i run over all space-time axes, sums over j run only over spatial axes, and 
sums over x run over all lattice sites. Above, A is a real-valued parameter irrelevant to the 
low-energy dynamics provided that the condition |A| > | is met. 



A very similar action appears to be first mentioned in [13], where Nielsen is credited for suggesting it. 
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In momentum space, we define — ^ ^^'^i^xi which renders the action in diagonal 
form, with the free Dirac operator given by 

( A;) = ^ sin - A ^ (cos kj-l). (2) 

There are precisely two zeros of this action for |A| > |, namely at the momenta k^^^ = 

(0, 0, 0, 0), and k^^^ = (0,0, 0, tt). Fields with restricted energy correspond to quark flavors. 
We choose the partition based on the ball,'^' B = {k^ : \k4^\ < |}, so that 

^V'^'^w- (3) 

On the other hand, careful consideration yields the second quark flavor^ 

^{k) =ia^^^^\T.Ak), (4) 

with shifted momentum given by Tj^^kf^ = (fc, A;4 + tt mod 27r). 

We can write the momentum-space action in terms of quark flavors using the partition 

in Eqs. (3) and (4). Deflning the isospinor ^(/c) = ( ^(2)^j^w' flavored action takes the 
form 

d'^k _ 



{2ny 



^{k). (5) 



^ i(7^ ® 1) sin kj^ - i\ (74 ® r^) ^(cos kj - 1) 
M j . 

Notice the momentum is restricted to B, and we employ a (spin ® flavor) notation rem- 
iniscent of staggered fermions. In addition to cubic invariance, the free action has the 
symmetries listed in Table I. The continuous vector and axial symmetries can alternately 
be described as a chiral symmetry: U{1)l^U{1)r. There is such chiral symmetry for the 
isosinglet, 1, and isovector, r^, flavor combinations. 

In the table, we have also spelled out how the coordinate space fleld, ip^, transforms under 
the momentum-space symmetries. The P, CT, T x flavor rotation, isovector vector, and 
isosinglet axial transformations are non-local in time. They depend on the time-smeared 
fields^ 

= C ^ ^ e*^"""' V.,.. , (6) 

and similarly for S-^ijjx, for which the energy integration is over the complement, namely 
B = {k^^: \k4\ ^ ^}. It is straightforward to verify that the coordinate space transformations 
in Table I leave the action in Eq. (1) invariant. Potential symmetry violating terms from 
the non-local transformations identically vanish because the momentum regions are non- 
overlapping, BnB — {}. 

* Throughout, we shall generically use energy to refer to the fourth component of momentum. 

^ The partition of the Brillouin zone into flavors is not unique, and alternate choices differ by lattice artifacts. 

Our partition is chosen in order to treat the flavors symmetrically. 
^ In the free theory, one is not required to define the second flavor in exactly this way. Such freedom, 

however, is lifted by the inclusion of gauge interactions, and Eq. (4) is the only possibility allowing for 

both fermion flavors to encounter the same gauge field. 
^ The temporal non-locality present in (5^ V'a; is of maximal range. A milder form of non-locality is needed 

in practice, and would result from a different definition of the flavors. For example, the infinitesimal 

transformation, ij^x — > V'a; + 1^75 (^x -1-4 + V'x-i) + • • • i corresponds to a isosinglet axial transformation 

different than that in Table I. While this transformation does not exponentiate simply, it is minimally 

non-local in time. 4 



TABLE I. Symmetries of the free minimally doubled action. Listed are the discrete and continuous 
symmetries of the action, and the corresponding transformations of the momentum-space and 
coordinate-space fermion fields. The spinor representations for the discrete symmetries are: C = 
7274; P = 74; and T = 7475, while the three-vector 1 = (1, 1, 1). Site reflection is an antilinear 
transformation, whereas the remaining transformations are all linear. 



Free Theory Symmetry 


Momentum Space 


Coordinate Space 






V'x- > 


Site Reflection 


e-'^^''^{T(g}l)W{k,-k4) 




P 


{P®l)^{-k,k4) 




CT 


{CT(^l)W{k,-k4) 


CT{5Bi^l-x4-h'^l-x4) 


T X Flavor Rotation 


{T (^icosOt'^ + i sin 9 r^)*(fc, -^4) 


ie--4 [e-'^5B^x,-x,-e'%i^^,-x,) 


U{1) Isosinglet Vector 






U{1) Isovector Axial 




e'^-'^^x 


U{1) Isovector Vector 




e'^^SB^x + e-^^i^x 


U{1) Isosinglet Axial 




e''^^5B^Px + e-''^^6^i^x 



B. Vector and Axial- Vector Currents 



For a generic current, Jij,{x), the lattice conservation law is precisely the statement: 
S//^^^Ai(^) — 0) where V* is the backwards difference operator, V* = — fx-ii- 
To derive the conserved currents, we appeal to the Noether procedure, which yields four 
conserved point-split currents. The two local symmetry transformations yield the isosinglet 
vector current 



(7) 



and the isovector axial current, J^^{x), which is given by the same expression as J^{x), but 
with the replacement: 7^ — )■ 7^75. The remaining two currents correspond to the non-local 
symmetry transformations in Table I. Applying the Noether procedure and appealing to 
the equations of motion yields the isovector vector current 



■^m(^) ^ \ ^B^^xlfi^B^x+fi + Ss^x+tiln^Bi^x - SQ^^-ff,SQ^a:+f, - % V'a; 

-ixs^j (^SBil^^-f^SBipx+j - SBi'^+fiiSBipx - s-^tp^-f^SsiPx+j + hi^x+juh'^x 



(8) 



and the isosinglet axial current, J^^^x), which is given by an expression identical to J^{x), 
but with the replacement: 7^ Ifilb- 



TABLE II. Symmetries of the gauged, minimally doubled action. The coordinate-space transfor- 
mations of the fermion and link fields are listed, along with the corresponding transformations of 
the free momentum-space fermion field. 



Gauge Theory Symmetry 



Coordinate Space 



Free Theory Momentum Space 




NaiVe CT 



Flavored Time Reversal (J-") 



Site Reflection (S) 
Naive Parity {V) 



^1-X,XA y-^j,l-X-J,X4^^4,l-X,X4, 



^x-X4, \y3,x-X4,-i f^4,a;,-a;4-l 



{CT®T^)W{k,-kA) 



C/(l) Isosinglet Vector 
U{\) Isovector Axial 



III. MINIMALLY DOUBLED FERMIONS WITH GAUGE INTERACTIONS 

A. Gauged action and its symmetries 

To include interactions with gauge fields, we gauge the isosinglet vector symmetry. On 
the lattice, the gauge fields are described by unitary link variables, C/^,x- In the presence of 
such a gauge background, the Wilczek action becomes 



The gauge theory retains only a subset of the symmetries possessed by the free theory. The 
non-local transformations listed in Table I no longer lead to symmetries because the link 
variables can carry momentum. Non-overlapping momentum regions are now connected 
by intermediate gluons, which can change quark flavor and chirality if the gluon energy 
is near tt. Symmetry transformations in the gauge theory are best described in coordinate 
space. Aside from cubic invariance and gauge invariance, the symmetries of the gauge theory 
are collected in Table II. While the continuous time reversal x flavor rotation symmetry 
is broken by gauge interactions, the particular flavor rotation hj 6 = ^ yields a discrete 
symmetry that we call flavored time reversal. Notice that the naive parity transformation, 
V, given by applying the parity operation to the lattice fleld, ip^, yields a symmetry of the 
interacting theory. This naive parity is not simply a parity transformation on each of the 
fermion flavors, but includes an additional isovector transformation by as well. A similar 
situation is encountered for the naive charge conjugation x time reversal transformation, CT. 
The transformation in Table II for the lattice fleld ipx, corresponds to a CT transformation 
in addition to an isovector rotation by r^. 

To arrive at an interacting theory with the correct continuum limit, we must augment 




(9) 
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the action^ with the following relevant and marginal terms [12, 15] 

>Cc = iCo *(74 (8) T^)* + Ci *'L>4(74 ® 1)* + C2 i^4ii^4i. (10) 

These Symanzik-level operators are allowed by the symmetries of the interacting theory. 
We assume the coefficient cq has been fine-tuned to remove the power divergence in the 
continuum limit, and that ci and C2 have also been tuned to recover SO (A) invariance. 

As the local, continuous transformations remain symmetries of the gauge theory, see 
Table 11, the corresponding symmetry currents are merely the gauged versions of Eq. (7). 
Applying the finite difference operator to the gauged symmetry currents, we find operator 
versions of the Ward-Takahashi identities. In an abritrary gauge background, the expecta- 
tion value of these divergences vanish, X]u(^u'^/i(^)) — Ylni'^t'^abi^)) — 0- 



B. Non-conserved currents 



The Noether procedure shows us the isovector vector and isosinglet axial currents have 
non- vanishing divergences, namely ^^<^^(^) = Dy{x), and V* J^5(a;) = Da{x), with 
the currents Jfi{x) and J,j,^{x) given by the gauged versions of Eq. (8), and their divergences, 
Dy{x) and Da{x), having the form 



(SS\U]\ (iS\U]\ 



\ S^x J B \ S^x Jb 

In the interacting theory, the temporal smearing is made gauge invariant by appending a 
Wilson line, yV4:{x,X4;x,y4), in the definition, Eq. (6). This Wilson line is the product of 
links in the time direction connecting the times and X4. Derivatives of gauged action 
appearing above have the form: (SS[U]/5iIj^) — D{x,x')5Bi>x', for example, where 
D{x,x') denotes the coordinate-space Dirac operator of the interacting theory. Notice we 
have the relations, Sefx + = fx, and Y.x9x{5Bfx) = Y.xi^B9x)fx-, for any functions 
and Qx- 

All operators in the divergences, Eqs. (11) and (12), connect different energy regions: B 
to B, and vice versa. Non-conservation of the isovector vector and isosinglet axial currents 
is hence due to flavor mixing. Without interactions, these divergences vanish due to flavor 
conservation in the non-interacting theory. Near the continuum limit, we can deduce the 
form of the anomalous divergences using the Symanzik effective action [16]. In this limit, the 
divergences, Dy{x) and Da{x), can be written in terms of continuum operators. The tower 
of such operators can be ordered by their dimension, with the lowest dimension operators 



In practice, of course, one adds the lattice versions of the terms in Eq. (10). For example, the term 
z^'(74 (8)r^)^ in the Symanzik action translates to the lattice operator i'tp^'y4ipx- The differences are of at 
least 0{a), and have no consequence on our computation of the chiral anomaly in the continuum limit. 
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TABLE III. Symanzik operator analysis of the anomalous divergences, Dy and Da- Transformation 
properties of the divergences, as well as the relevant and marginal operators in the Symanzik 
effective theory are listed. For simplicity, we have suppressed all coordinate dependence by listing 

the transformations as either even or odd. The coordinates of all operators transform in a covariant 
manner under the listed symmetries. We have listed operators that are covariant under cubic 
rotations, and invariant under U{l)v- Operators not invariant under the isovector U{1)a must 

be accompanied by an explicit power of the quark mass, ruq. The arrowed covariant derivative is 

^ <- 

defined hy D = D* — D, with Z)^ = + iA^. Notice we have utilized the equations of motion to 
eliminate ]p terms. 







Transformation 




Operator 




V 


CT 






— 


+ 


— 


— 


*(74 <8) 1)^' 


- 


+ 


+ 


— 


1^(74 (8)r2)^' 


- 


+ 


+ 


+ 




— 


+ 


+ 


+ 






+ 


+ 




-D4['*f74 (8> 




+ 




+ 


-D4[*(74<8)T^)*] 




+ 






imq^{l (g) 1)^' 




+ 


+ 


+ 






+ 


+ 








+ 


+ 


+ 






+ 


+ 


+ 


Da 










*(7475 ® 1)* 






+ 


+ 


*(7475 (8> T^)* 






+ 




z^'£)4(7475 (8) 1)* 






+ 




i*£)4(7475 (8>r^)^' 






+ 


+ 


Z?4["^(7475 ® 1)*] 










L>4[«'(7475 ®T''^)^] 








+ 


mq^{-f5 (g 1)^' 










"lg"^(75 (gr^)^ 








+ 













as the most relevant in the continuum limit. Our concern is only with those operators that 
survive when the lattice spacing, a, is taken to zero. 

In Table III, we list the transformation properties of the divergences, Dy and Da, under 
the discrete symmetries of the lattice theory. We also list all possible relevant and marginal 
operators classified as to their transformation properties. From the table, we see that the 
divergence of the isovector vector current is required to have the form 

D3 ^cyL'4[^(74(8T3)*] + C»(a). (13) 
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In other words, conservation of the isovcctor vector current in the continuum hmit requires 
the fine tuning of a marginal operator. As the marginal operator breaks 5*0(4) invariance, 
one might think that tuning the action to recover 5*0(4) invariance in the continuum, via 
the terms in Eq. (10), would be enough to render cy — 0. This, however, is not the case. 
While the divergence of the isovector vector current depends on functional derivatives of 
the renormalized action, Eq. (11), there is additional 5*0(4) breaking necessitated by the 
temporal non-locality of fiavors. Even with a more local definition of the isovector vector 
current, ultimately 5*0(4) must be broken, and fine tuning will be required so that Dy{x) — 
in the continuum. 

The divergence of the isosinglet axial current is required by discrete symmetries to have 
the form 

Da = CAi mg^{-f5 ® 1)* + CA2 ef,uai3Ff,^F^0 + cas L>4[*(7475 ® 1)*]. (14) 

The coefficient of the first term, is fixed by operator algebra, Cai = 2, while the coefficient 
of the second term, ca2, should be fixed by the chiral anomaly. We pursue the computation 
of the anomaly in the next section. The last term is a surface term not excluded by the 
lattice symmetries, and modifies the divergence of the axial current in the chiral limit. This 
5*0 (4)-breaking operator is similar to that encountered in the divergence of the isovector 
vector current, Eq. (13). Even with an action fine tuned to recover 5*0(4) in the continuum, 
the coefficient ca3 will not in general vanish. This is due to 5*0(4) breaking introduced 
by the temporal non-locality required to form the isosinglet fiavor combination in the axial 
current. Alternate definitions of the current must share this malady, and one must fine tune 
the current to arrive at ca3 = 0. We will assume this tuning has been performed, so that 
the operator e^^a^F^uFap (or e^i^F^y in two dimensions) is the only contribution to Da in 
the chiral and continuum limits. 



IV. CHIRAL ANOMALY 

For perturbatively weak gauge backgrounds, we can compute the divergence of the isosin- 
glet axial current using lattice perturbation theory. It is useful to employ a compact notation; 
and, to this end, we write the minimally doubled action in the form, 5* = ^ il)^{x\D\y)il)y^ 
so that, as an abstract operator, the propagator is merely D~^. Prom Eq. (12), the expec- 
tation value of the axial current's divergence in a gauge background is 

{Da{x)) = Tr75(a; 

where we have used {75, D} = 0, and the trace is over spin. The flavor-changing propagators 
are given by Wick contractions of time-smeared flelds 

D^l{x,x') = {6-gilj^6B^^,), (16) 

with a similar expression for D~^{x,x'). Expanding in powers of the external gauge field, 
we write D = D^'^^ + D^^^ + . . ., where the ellipsis represents terms with at least two powers 
of the gauge field. 

For simplicity, we restrict our attention to the Abelian anomaly in two dimensions. A 
further simplification is afforded by adopting Coulomb gauge, for which the time-links are 



\x) - (B^B), 



(15) 
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FIG. 1. Leading contributions to the two-dimensional Abelian anomaly in lattice perturbation 
theory. The squares denote insertion of the isosinglet axial current's divergence, while the circle 
denotes the gauge vertex. Fermions are directed lines, gauge fields are curly lines, and the flow of 
momentum from the axial current is shown with dashed lines. 



unity. The leading contributions to the divergence appear at linear order in the gauge field. 
Thus we have 



(17) 



The first trace corresponds the first diagram in Fig. 1, while the second trace corresponds 
to the second diagram. Both diagrams, however, are identical, and are easiest to evaluate 
in momentum space. 

In momentum space, the gauge vertex, M.{k -\- q,k), in Coulomb gauge is defined by 



(18) 



Using the Wilczek action, we have 

The free Wilczek propagator in momentum space has the form 

dW(A;) = -L»W(A;) V-^{k), 



(19) 



(20) 



with D^^\k) = i'jismki + i72[sinA;2 — A(cos /ci — 1)], which is the two-dimensional version 
of Eq. (2), and T>{k) = sin^ /ci + [sin/c2 — A(cos/ci — 1)]^. Injecting momentum q into the 
divergence, we have the momentum-space result 



{DA{q)) = 2 



d^k 



rTr75 



D^^^~\k) + {k + q)\M{k + q,k)-{B^ B).(21) 



keB, k+qeB 



Contributions to Eq. (21) around = (0,0) appear only at second order in q. The 
minimally doubled action, however, possesses two low-energy regions, and there are linear 
contributions around the momentum = (0,7r). Writing q^ — )> 116^2 + q'/j,, we then expand 
q'^ about zero retaining all linear terms. First, we notice that restrictions on the integration 
region allow us to split the integration into a finite region and an infinitesimally small region. 



I -I- I ■ 

keB,k+q'eB keB keB,k+q'eB 



(22) 
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The integral over the infinitesimal region can be shown to vanish quadratically in q'^. Second, 
we notice that in the remaining finite interval, the integrand is odd in k2 if ^2 — 0- Thus 
the only linear contribution is proportional to q'21 and, after some algebra has the form 

^Tr75 ^^^^^ M{Kk). (23) 

A valuable simplification is shifting the energy integration by tt in the {B ^ B) term. At 
this stage, we utilize Tr 757^71, = 2ie^^ to perform the spin trace. The divergence of the 
axial current then takes the form 

s./n/" cPk d I A sin^ /ci — cos fcifsin A;2 — A(cos fci — 1)1 



(Pk d / sin A';i COS A:2 \ , 
= _«efe-.MiM/p^g^(^i^). (24) 

jfcGB 

As the surface terms vanish, the only contribution to the divergence arises from the region 
where \ki\, |/c2| < £, with £ — >■ 0. The integrand has an integrable singularity in this region. 
Carrying out the integration for finite e, and then taking the limit £ — > produces the 
correct chiral anomaly 

2ie eNf 
{DA{q)) = (^2 - 7T)A,{q) --^e F^,, with Nf = 2. (25) 

V. MINIMALLY DOUBLED WEYL FERMIONS 
A. Free Action 

Using the chiral projector, Pr = |(1 + 75), on a minimally doubled action, we can reduce 
the number of spin components by half. The resulting action can be written in terms of a 
two-component Weyl spinor Xx, namely 

^^\Y1 [xU'^uXx+ix - xl+nicr^Xx - y J] [xlXx+j + xl+jXx - '^xlxx] , (26) 

X,IJ. x,j 

with = (cr, —i). There are again two low-energy modes, which we now label hj kli — 
(0,0,0,0) and ki^^ = (0,0,0,7r). Fields with chirality can be created by partitioning the 
Brilluoin zone. We take xik) — x'^^Hk), for the right-handed field, and x(^) = 

X^^^T-jj-^k), for the left-handed field. Forming a Dirac spinor from the left- and right-handed 

fields, ip{k) = and its conjugate spinor, ijj{k) — — X^^^^(A:), X^*-^^(/i;)^ , we can 

write the action in the form 

f d^k — [ \ ^ 1 

S ^ j^j^^'4^{k)Y^i-^,,smk^-i\'-^Al^^{coskj -I) ip{k), (27) 
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with Euclidean chiral basis Dirac matrices 7^ 







-la 



ia^ 



{a,i), and 75 



diag(l,-l). 

The free momentum-space action has: cubic symmetry; invariance under the product of 
charge conjugation and parity, CP; time reversal invariance, T; U{1) vector symmetry, e*^; 
U{1) axial symmetry, e*^'*'". In coordinate space, the continuous vector symmetry is local, 
while the axial symmetry is non-local in time 



Xx 



U{l)y 

e^%X. + e-%X., U{1)a 



(28) 



Thus the lattice current corresponding to the U{l)v symmetry. 



Xx'^'^nXx+fj. ''r Xx+fj.'^^IJ.Xx '^X^fJ.jiXxXx+j Xx+jX^ 



\ A.xA.x+3 y^x+j 

is conserved, while that corresponding to the non-local U{1)a symmetry. 



(29) 



■^^5(2;) = - Ssxlif^^i^BXx+i, + SBxl+^,^(^^^SBXx - % xl;^t^M% Xa;+/. - S^xi+^,^(^^lSBXx 

-iXS^j (^Sb xl^B Xx+j - Sb xI+j^b Xx - h xlh Xx+j + % xl+jh X^ , (30) 
is only conserved in the free theory. 



B. Gauge Theory and Chiral Anomaly 

The inclusion of gauge fields is straightforward; one gauges the U{l)v symmetry of the 
free theory. As in the case of the Wilczck action, additional relevant and marginal terms 
must be added with fine-tuned coefficients in order to arrive at the correct continuum limit. 
As the terms needed are similar to those given above in Eq. (10), we do not spell them out 
explicitly. 

The vector symmetry transformation in Eq. (28) acts locally on the fermion field, it re- 
mains a symmetry of the interacting theory. The corresponding vector current is the gauged 
version of Eq. (29), and has vanishing divergence in an arbitrary gauge field background. 
The same fate does not befall the axial symmetry in the presence of gauge interactions. The 
gauged axial current in Eq. (30) has non-vanishing divergence 

i:v;...^(fa)^..-(^)^^..^..(5)^-..Km- 

(31) 

Left-handed and right-handed fields mix due to interactions with gauge fields, and this 
mixing is the source of the chiral anomaly. 

Due to the breaking of SO {4) invariance that is inherent to defining the axial current, 
Eq. (30), fine tuning is required to recover the chiral anomaly in the continuum limit. At 
the Symanzik level, the divergence of the axial current receives a contribution from an 
5*0(4) breaking operator, namely D^lipj^j^ip]. We assume that this contribution has been 
tuned to vanish by improving the current. This leaves only a potential contribution to 
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the divergence from the pure gauge operator e^^vapF^vFap (or, in two dimensions. e^„F^i,). 
Using the divergence of the axial current, Eq. (31), we can compute its expectation value 
in a gauge field background to determine the coefficient of the pure gauge operator. For 
simplicity, we work in two dimensions, for which the Weyl spinors are one-component objects, 
and cr^ = (1, —i). In a weak gauge field, we arrive at the diagrams shown in Figure 1. Their 
evaluation parallels that above in Sec. IV, and one encounters integrals identical to those 
arising in the computation of the axial anomaly, Eq. (24) . Notably absent is the trace over 
spin indices, which reduces the overall value by a factor of two. Consequently we find. 



Thus the continuum limit of the minimally doubled Weyl action, Eq. (26), describes a single 
Dirac fermion, and reproduces the chiral anomaly. Fine tuning of both the action and the 
axial current, however, is required to obtain this result. 

VI. SUMMARY 

Through the investigation of minimally doubled fermions, we have explained the ap- 
pearance of the axial anomaly in chirally symmetric lattice actions. Although such free 

actions have a fiavor-singlet axial symmetry, the lattice field transforms non-locally under 
this symmetry. When gauge interactions arc included, the symmetry is broken by fiavor 
changing mediated by gluons carrying the requisite momentum. In the continuum limit, we 
demonstrated that these interactions produce the correct chiral anomaly in two dimensions. 
Similar fiavor-changing interactions lead to non- vanishing divergences of non-singlet vector 
currents. To arrive at a conserved non-singlet vector current, fine tuning is required to 
restore SO {A) invariance. The singlet axial current must also be tuned to recover 5*0(4). 
Finally we considered the chiral projection of minimally doubled fermion actions. Without 
gauge fields, these projected actions describe a single, chirally symmetric Dirac fermion. 
Chirality-changing interactions are present in the gauged theory; they supply precisely the 
right contribution to the chiral anomaly. 
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WILSON LATTICE ACTION 

It is straightforward to compute the chiral anomaly for the Wilson action in momentum 
space. We present a simple derivation. The Wilson action takes the form 




(32) 



Sw[U] = 2 SI \^xU^^,x{.lu^ + r)il^x+„ - ip^+^UlAlf^ - r)ip^ - 2(r - mo) i'^ipx , (33) 
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where r and mo are real- valued parameters, the latter is the bare value of the fermion mass. 
Because of chiral symmetry breaking, the bare mass receives an additive renormalization 
proportional to as /a, where a is the lattice spacing. Below, we implicitly tune the bare 
parameter mo to exactly cancel the power divergent mass renormalization. This fine tuning 
is required to arrive at the correct anomaly in the continuum limit. 

Even with vanishing renormalized quark mass, the axial transformation, ip^ — )■ e'^^'^^ipx, 
is not a symmetry of the Wilson action. As a consequence, the axial- vector current, J/^six), 
given by 

1 



(34) 



is not conserved. Writing the action as Sw — y^x{^\^w\y)i^y, with the Wilson-Dirac 
operator defined by Dw — Ip-\- R, casts the divergence in the form 



(yiJM) = -Tr-i^{x\[D^,R]\x)=T:Y-f^{x\[D^\lp]\x) 



(35) 



Expanding in powers of the gauge coupling, we write = D^^"* + Z)^^^ -|- . . ., and R — 
+ -R*^^-* + . . .. In two dimensions, the leading contribution to the divergence occurs at 
first order in the gauge field 

{KW)-^i^{d'T\iP''']-{d'T\P''' + R^^^^^^ (36) 

The first anti-commutator is graphically depicted as the first diagram in Fig. 1, while the 
second anti-commutator is depicted by the second diagram. In momentum space, the free 
Wilson propagator is 

D^w~\k) = J]H7^sinA;^ + r(cosfc^ - l)]V-\k), 



(37) 



with T){k) = J2u si^^ + [S;/(cosA;j/ — 1)]^, and the gauge interaction vertices are 
5R^^^/SA^ ^n^{k + q, k) = ^[e*'^'^ - e'^'^'^+i^^ 



Evaluating the traces, we find 



(Pk 



. (27r)2 



dk. 



cos kfj, sin k^ 



+ 



d ( cos ky sin k,, 



dk,, 



V{k) 



(38) 
(39) 

(40) 



where wc have kept the two terms ordered as to their corresponding diagrams. Non- vanishing 
contributions to the integral arise from the region, \kv\ < e, where the integrand has 
an integrable singularity. Taking £ — > 0, we obtain 



eNf 

-^e^,F^,{q), with Nf = 1. 



(41) 
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